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PROPAGATION  EFFECTS  OF  EMC-RESISTANT  SYSTEMS 
IN  COMMUNICATION  AND  NAVIGATION 


THE  ATMOSPHERIC  PROPAGATION  MEDIUM  BETWEEN  45  AND  75  GHz 


Hans  J.  Llebe 

Institute  for  Teleconmunication  Sciences 
National  Telecommunication  and  Information  Administration 
U.S.  Department  of  Commerce 
Boulder,  Colorado  80303 

SUMMARY 

The  clear  atmosohere  (dry  air  mass)  represents  a  unique  filter  over  the  45  to  75  GHz  range  with  fre¬ 
quency-dependent  absorption  properties  caused  by  the  microwave  spectrum  of  oxygen  not  found  at  any  lower 
frequency.  The  basic  information  necessary  for  determining  molecular  line,  band,  and  continuum  absorotion 
is  contained  in  a  transfer  function  model.  Close  to  sea  level,  an  unstructured  band  exists,  centered  at 
60  GHz  with  a  half-width  of  8  GHz  and  a  pressure-dependent  intensity  that  is  corrected  for  pressure- 
induced  line  overlap  effects.  At  elevations  above  12  km,  the  lines  separate  and  radio  channels  with  up 
to  400  MHz  bandwidth  can  be  accommodated  between  lines.  Above  30  km,  the  line  properties  become  polar¬ 
ization  and  orientation  (i.e.,  with  respect  to  the  direction  of  a  magnetic  field)  dependent  due  to  the 
Zeeman  effect.  Each  oxygen  line  splits  proportionally  with  geomagnetic  field  strength  (0.2-0. 3  Gauss) 
into  numerous  sub-lines,  wmch  conform  to  Zeeman  patterns  extending  over  a  megahertz  scale  around  eacr 
1 i ne  center. 

Using  this  model  m  conouter  routines  in  conjunction  with  distributions  of  pressure,  temperature,  a no 
humid! ty  along  a  radio  patn  makes  it  possible  to  predict  frequency-dependent  signal  attenuation  and  delay 
effects  of  the  neutral  atmosphere  (h=0  to  ICO  km).  Further  analysis  investigates  possible  Dandwidth 
limitations  of  the  propagation  medium  by  employing  short-pulsed  (<_  1  ns)  test  signals  with  carriers  in  the 
45  to  75  GHz  range.  Examples  are  presented  of  specific  attenuation  (dB/km)  and  excess  delay  (ps/km)  rates, 
of  zenith  path  attenuation  (dB),  and  of  pulse  distortion  over  a  terrestrial  link  (50  GHz,  35  km)  for  Gaussian, 
rectangular,  and  BPSK  modulation.  Benefits  from  the  frequency-variable  atmospheric  properties  include  a 
fairly  well  defined  signal  range,  affording  transmission  security  and  trade-offs  between  range  and 
susceptibility  to  detection  and  interference. 


1.  INTRODUCTION 

The  clear,  nonturbulent  atmosphere  sets  limits  on  the  transmission  of  radio  waves  in  the  45  to  75  GHz 
frequency  interval.  Absorption  due  to  0?  causes  signal  attenuation,  changes  in  transit  time,  and  increases 
in  the  background  noise.  Accurate  knowledoe  of  the  atmospheric  transfer  function  provides  the  basis  for  a 
variety  of  specialized  applications  including: 

.  Broadband  communication  with  a  restricted  and/or  frequency-dependent  radio  path  range; 

.  Shielding  of  satellite-to-satellite  links  from  ground  interferences; 

.  Secure  communication  between  high-flying  aircraft  and  alrcraf t-to-satel 1 i te  links; 

.  Remote  sensing  of  vertical  temperature  profiles  from  satellite  or  ground-based  observation 
platforms. 

In  the  past,  the  accuracy  of  predicting  atmospheric  transfer  properties  from  meteorological  data  was  poor. 

ITS  developed  a  new  spectrometer  which  measures  dispersion-pressure-profiles.  These  were  reduced  to  reliable 
spectroscopic  coefficients  for  the  more  than  30  oxygen  absorption  lines  centered  around  60  GHz  [LIEBE  et 
al.,  1977],  The  extensive  laboratory  data  reduce  errors  in  predicting  attenuation  rates  over  the  range 
0.1  to  20  dB/km  to  a  few  percent  (1  to  3").  The  standard  microwave  refractivity  with  two  coefficients 
(dry  and  wet  term)  was  expanded  Into  a  250-coefflcients-scheme  computer  model  [LIEBE,  1981],  This  paper 
discusses  the  modeling  strategy  and  provides  a  series  of  mostly  graphical  outputs,  which  elucidate  the 
influence  that  frequency,  range,  humidity,  altitude,  and  geomagnetic  field  strength  have  upon  transfer 
properties  of  the  neutral  atmosphere  over  the  height  range  h=0  to  100  km. 


2.  ATMOSPHERIC  TRANSFER  PROPERTIES 

2.1  Transfer  Function  and  Specific  Propagation  Rates 

Amplitude  and  phase  response  of  a  planar  radio  wave  traveling  the  distance  L  in  km  and  starting  with  the 
field  strength  EQ  are  described  by  the  transfer  function 

T  *  (E/Eq)  *  exp  ( TL)  «  exp  [j0.02096f(106  ♦  N)L]  (1) 

where  frequency  f  is  In  gigahertz  (GHz)  and  r  In  km"^  is  the  propagation  constant. 

Germane  to  any  propagation  model  is  a  convenient  macroscopic  measure  of  the  interaction  between  radiation 
and  the  atmosphere.  Complex  refractivity  N  (in  parts  per  million),  expressed  in  terms  of  measurable 
quantities,  assumes  that  role.  For  air,  N  consists  of  three  components 

N  *  Nq  ♦  0(f)  ♦  JN" (f )  ppm  (2) 

namely,  the  frequency  independent  refractivity  N0  plus  various  spectra  of  refractive  phase  dispersion 
0(f)  and  absorption  N"(f).  Usually,  the  imaginary  part  of  Eq.  2  is  expressed  as  the  specific  power 
attenuation  a  and  the  real  part  determines  the  phase  delay  U  (with  reference  to  vacuum)  as  follows: 
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a  =  0. 1820f  N" 


d8/km 


and 


8  =  0.02096f(NQ  +  0)  radians/km 


Accordingly,  the  propagation  constant  r  and  the  excess  propagation  delay  time  t  are 
r  =  -0.1151a  +  j ( 2 . 096  - 1 04f  +  d)  1/km 
t  *  (d/2rf )  103  =■  3.336(Nq  +  D)  ps/km. 


(3) 

(4) 

(5) 

(6) 


2.2  Complex  EHF  Refractlvlty  of  Air 

The  frequency-independent  refractlvlty  Is  given  by  [HEBE,  1981] 

Nq  *  (2.589p  +  41.6e8  +  2.39e)0  ppm 

and  can  be  calculated  in  a  straightforward  manner.  A  physical  state  of  air,  which  determines  the 
absorber  population  in  a  radio  path.  Is  described  by 

p  dry  air  pressure  (barometric  pressure  P  =  p  +  e),  kPa  (1  kPa  =  10  mbar); 

0  *  300/T  relative  inverse  temperature  (T  in  K); 
e  water  vapor  partial  pressure,  kPa;  and 

w  cloud  water  concentration,  g/m2 


(7) 


(8) 


The  absorption  and  dispersion  spectra  are  formulated  from  line  contributions,  a  continuum  N  ,  and  a  cloud 
water  term  Nw;  i.e. , 

N"(f)  =  £  (SFM)i  +  N“c  +  N“w  ppm  (9) 

1 

and 

o(f)  -  £  (SF')i  +  N'c  +  N'w  ppm  (10) 

Both  expressions  are  discussed  in  detail  by  LIEBE  [1931]. 

Oxygen  absorption,  which  dominates  in  the  45  to  75  GHz  range,  is  calculated  from  the  line  strength 


S  =  a^d3  exp  [a2(l  -  9)] 
and  the  shape  factors  F"(f)  and  F'(f). 


kHz 


(ID 


For  the  line  shape,  we  chose  the  VanVleck-Weisskopf  function  as  modified  by  ROSENKRANZ  [1975]: 


F"’^[ 


y  -  (v0  -  f)6 
(v0  -  f)2  ”  y2 


y  -  (vQ  +  f)6  -I 

(v0  +  f)2  *  yZ  J 


A  Hilbert  transform  of  (12)  yields  the  dispersion  profile  [BROUN,  1982] 


GHz' 


(12) 


F' 


(y0  -  f)  *  y  (y  +  8f)/v0  (vQ  +  f)  +  y(y  -  sf )/vQ 


<v0  -  i)2  +  y2 


(v0  +  f)2  ♦  yZ 


GHZ 


-1 


which  is  slightly  different  from  F'  given  by  LIEBE  [1981],  The  line  parameters  are 

.0.9 


Y  *  a3(p  +  1.3e)e 
&  *  a4pea5 


GHz 

1 


(13) 


(14) 

(15) 


where  y  is  the  width  and  |6|  <  1  Is  the  overlap  Interference.  The  spectroscopic  coefficients  a,  .  and 
the  molecular  line  center  frequencies  v  for  1  =  46  0.  lines,  as  well  as  for  i  »  30  H.O  lines,  ° 
are  tabulated  by  LIEBE  [1981].  0 

The  continuum  spectrum  consists  of  a  water  vapor  and  a  dry  air  absorption  term, 

N"c  ’  fpe3(19e  ♦  74p/[f2  ♦  (0.012p9)2])10'6  ppm  (16) 

while  refractive  continuum  phase  dispersion  Is  neglected,  that  is 

N*c  *  0  (17) 

Liquid  water  extinction  N”v  of  suspended  droplets  (cloud,  fog,  haze,  aerosol  hydrometeors)  can  be  derived 
from  dielectric  data  of  bulk  water  given  by 


* - * . .  - . 


and 


c'  »  4.9  *  (185.1  -  1 1 3/9 )/[l  +  (fr)2] 


and 

c"  =  ft(185. 1  -  113/e)/(l  ♦  (fx)2] 

;J3) 

where 

r  =  0.4170  10‘4  exp  (7.136)  ns 

(19) 

The  Rayleigh  absorption  approximation 

N"w  =  4.49w  £”/[(£■  +  2)2  *■  U")2] 

{20) 

of  Mie  scattering  losses  is  appropriate,  and  refractive  dispersion  is  neglected  (i.e.,  N'  *  0). 

w 

2.3  Radio  Patf  Transfer  Function 

While  equations  (3)  and  (4)  describe  the  specific  rates  i  and  3  of  a  homoqeneous  air  volume  element  over  a 
path  increment  dx,  the  propagation  medium  of  an  extended  radio  path  generally  is  inhomogeneous.  Total 
path  attenuation  A  and  excess  delay  q  over  a  length  L  are  evaluated  as  follows: 

L 

A  »  /  a(x)dx  d8 

o 

and 

,  (2D 

0=7  S(x)dx  radians 

o 

Equation  (21)  allows  us  to  define  a  radio  path  transfer  function 

Tl  =  exp  (-0.115A  +  jo)  (22) 

The  excess  propagation  delay  time  of  the  path  Is 

tL  =  0. 1 59t/f  ns  .  (22a) 


3.  ATMOSPHERIC  TRANSFER  PROPERTIES  FOR  ALTITUDES  0  to  30  KM 

The  complex  refractlvlty  N  (Eq.  2),  as  formulated  In  Section  2.2,  is  evaluated  by  a  computer  routine,  wmch 
produces  values  of  attenuation  a(f),  dispersion  D(f),  and  refractivity  N0.  The  program  is  valid  for  fre¬ 
quencies  up  to  1000  GHz  and  for  heights  up  to  30  km.  In  the  troposphere  (h  <  10  km),  only  lines  of  H,0  and 
O2  are  important.  The  calculation  of  fl(f,  P,  T,  e,  w)  is  done  in  a  subroutine  usino  line-by-line  super¬ 
position  at  each  frequency,  the  continuum  (Eq.  16),  and  the  suspended  hydrometeor  absorption  (Eq.  20). 

Water  vapor  pressure,  e(kPa),  absolute  humidity,  v  (g/m3),  or  relative  humidity,  RH(t),  can  be  specified. 

The  main  program  generates  a  frequency  grid  with  a  resolution  that  is  tjed  to  the  line  center  frequencies 
which  fall  within  the  des'^J  range.  Two  consecutive  values,  Vq  and  ug,  are  taken  and  the  number  of  points 
between  them  is  specifies  -.jally  2  to  5).  Also  automatically  added  to  the  grid  are  the  half-power  points 
at  vo  -  Y-  Thus  it  is  a  1  that  attenuation  maxima  and  dispersion  peaks  will  be  plotted.  The 

program  feeds  one  freque-  :  a  time  into  the  subroutine  and  receives  values  of  a  and  D,  which  are 

stored  temporarily  in  two  es.  When  all  numbers  of  the  grid  have  been  processed,  the  output  is  presented 
in  a  numerical  printout  and  orocessed  by  a  graphics  routine.  Examples  are  depicted  in  Fiqs.  1  to  3 
and  Table  1,  taking  w  *  0. 

An  overview  of  the  full  program  to  1000  GHz  at  sea  level  is  given  in  Fig.  1.  Significant  attenuation 
extends  over  several  orders  of  magnitude.and  above  120  GHz  is  caused  mainly  by  water  vapor.  The  low- 
frequency  wing  of  the  dispersion  spectrum  from  the  strong  557-GHz  H20  line  makes  its  influence  known  down 
to  about  50  GHz.  The  frequency  range  45  to  75  GHz  is  dominated  by  oxygen  absorption,  as  evident  from 
Figs.  2  and  3.  The  change  In  attenuation  from  h=0  to  6  km  is  small.  Above  hjlO  km,  however,  one  can  see 
a  significant  difference  -  the  unstructured  60  GHz  band  shape  transforms  into  a  highly-structured  single 
Tine  response  with  distinct  peaks  and  valleys.  The  lines  narrow  with  increasing  altitude,  deepening  the 
valleys  while  the  line  peaks  stay  approximately  constant.  This  behavior  is  illustrated  in  Fig.  3  in  more 
detail  and  backed  by  selected  numerical  examples  in  Table  I.  Tracing  the  line  peaks  to  zero  intensity 
over  the  range  h»30  to  100  km  introduces  complications  of  Zeeman  splitting  by  the  weak  geomagnetic  field 
(0.3  to  0.6  Gauss),  wh^-h  Is  discussed  In  the  following  section. 

Frequency-agile  systems  which  operate  within  the  60  GHz  band  can  have  the  capability  of  selecting  a 
desired  atmospheric  path  attenuation  to  limit  range  or  afford  interference  protection.  Fixed  rates  up  to 
lJ-  '°  dB/km  can  be  realized  below  h*8  km,  and  they  drop  to  .  <  1  dO/km  at  h=16  km.  Another  property 
of  the  60  GHz  band  is  a  frequency-dependent  temperature  dependence  of  ,  and  Applying  the  approximations 

a(T)  =  a(300K)ex  and  S(T)=  8(300K)9y  (231 

to  sea  level  air,  we  obtain  [UEBE  and  HOPPONEN,  1977]  (a, 8  -  see  Table  m  and  Fig.  2): 
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f 

50 

55 

60 

65 

70  GHz 

X 

2.49 

1.30 

2.35 

1.20 

2.71 

y 

1.84 

1.98 

3.30 

1.87 

1.88 

4.  MESOSPHERIC  PROPAGATION  MODEL 

The  model  calculates  specific  attenuation  of  a  single  line  above  h  =  30  km.  Decreasing  pressure  in  Eq.  (14) 
causes  it  to  eventually  approach  the  finite  Doppler  line  width 

Yq  =  6.20vq/  ^  md  kHz  (24) 

where  m  is  the  molecular  weight  (e.g.,  O2.  m  =  32,  ‘3  =  1.  v  =  60  GHz;  =  65.8  kHz).  An  adequate  approxi¬ 
mation  can  be  made  by  replacing  y  with 

Yh  =  >'  *  ♦  (25) 

Zeeman  splitting  of  isolated  oxygen  lines  due  to  the  influence  of  the  earth's  magnetic  field  strength  H 
introduces  considerable  complications.  Each  line  is  identified  by  a  quantum  number  K.  The  presence  of 
the  steady  field  H  splits  K1  lines  Into  three  groups  of  ( 2K-  +  1)  sublines,  thereby  redistributing  line 
attenuation  over  a  fixed  frequency  range.  In  principle,  three  different  Zeeman  patterns,  a.  ,  ,,  are 
possible  for  any  K1  line.  These  patterns  are  obtained  by  evaluating  n  and  :  coefficients  in’"’ 
it  and  n-  groups  as  detailed  by  LIEBE  [1981],  Each  sum  of  sublines  determines  a  Zeeman  pattern 
according  to 


“1,2,3  =  “0  &  {1  +  [(f  *  vo^h]2} 
where  v*  Is  a  function  of  n  and 

a0  ■  O.182vo(a1p/Yh)03exp  [a2 ( 1  -  6)]  dB/km 


(26) 


(27) 


Additional  information  has  to  be  provided  before  the  mesospheric  O2  line  pattern  tz  can  be  calculated,  that 
is,  the  polarization  of  the  radio  wave  and  its  orientation  with  respect  to  the  neomannetic  field  H.  For 
linearly  polarized  radiation,  an  angle  of  orientation  O  is  defined  between  the  wave's  magnetic  field  com¬ 
ponent  in  a  plane  of  constant  phase  (i.e.,  perpendicular  to  the  direction  of  propagation)  and  the  geomagnetic 
field  direction  H.  Line  attenuation  then  generally  consists  of  a  mixture  of  the  three  Zeeman  patterns: 

7  2 

az  *  0^  cos  0  +  (a2  +  a-j)  sin  0  (28) 

Mesospheric  Oj  line  properties  are  direction  (anisotropic  medium)  dependent,  as  implied  by  Eq.  (28). 

Further  coordinate  transformation  is  necessary  when  a  fixed  antenna  receives  radiation  from  a  mesospheric 
space  element  [HARTMANN,  1981], 


The  patterns  a.,,,,  of  the  O2  lines  K  =  1"  to  29'  have  been  calculated  for  local  p-T-H  conditions  over  the 
height  range  30  to  100  km  at  two  magnetic  field  strengths,  H  =  0.3  and  0.6  G.  Two  examples,  K  =  7+  and  19', 
are  depicted  In  Fig.  4.  The  complete  set  K1  «  1  to  29  is  available  upon  request  (to  be  published  in  a  report). 
Above  h  *  70  km,  the  individual,  now  mostly  Doppler-broadened  components  become  di^cernable.  Although  details 
depend  on  specific  K  values,  the  relative  features  at  given  heights  are  similar  to  the  cases  illustrated 
in  Fig.  4. 


The  results.  In  combination  with  beginning  dissociation  to  atomic  oxygen  (e.g.,  at  h  =  120  km,  0?/0  *  0.25), 
establish  h  *  100  km  as  the  plausible  boundary  to  outer  space  for  radio  path  modeling.  The  cutoff  height 
for  water  vapor  is  much  lower  (h  «  20  km). 


5.  ZENITH  PATH  ATTENUATION 

Cumulative  transfer  characteristics  of  the  neutral  atmosphere  are  evaluated  by  a  numerical  integration  of 
Eq.  (21).  For  a  zenith  path,  the  increments  dx  are  equal  to  dh,  and  dh  is  chosen  so  that  each  slab  is 
quasi -homogeneous  (p,  T,  e,  and  w  are  constant).  The  formalism  developed  in  Sections  2  to  4  is  applied  to 
calculate  a(dh)  and  8(dh).  Height  profiles  of  the  meteorological  quantities  are  obtained  from  model 
atmospheres  or  in  situ  data  (e.g.,  radio-sonde).  Total  precipitable  water  vapor  and  hydrometeor  absorption, 
plus  the  relatively  stable  dry  air  attenuation,  are  the  determining  quantifies  jf  the  radio  path  transfer 
function  (Eq.  22).  The  water  vapor  and  cloud  water  column  densities)!  g/m-’km  =  1  mm)  are  accounted  for  by 

20  km  10  km 

V  *  Z.219  /  e(h)0(h)  dh  mm  and  W  *  f  w(h)  dh  ntn  (29) 

0  0 

Ory  air  attenuation  Is  evaluated  for  all  O2  lines  up  to  h  *  30  km  (Sections  2  and  3),  and  for  isolated  O2 
lines  over  the  frequency  range  v0  +  100  MHz  by  considering  the  Zeeman  effect  (Section  4). 

A  standard  example  is  the  one-way  zenith  response  through  the  U.S.  Standard  Atmosphere  [NOAA,  1976], 

Zenith  attenuation  is  displayed  in  Fig.  5,  where  the  calculation  was  performed  "layer-by-layer",  encompassirg 
53  slabs  between  h  >  0  and  100  km.  One  notices  that  outside  the  O2  absorption  frequency  ranges  (50  to  70  GHz 
and  115  to  123  GHz),  water  vapor  V  plays  the  dominant  role.  The  frequency  range  53-  to  67  GHz  is  detailed 
In  Fig.  6  assuming  conditions  (H  *  0.3  Gauss,  m-pattern)  which  yield  the  maximum  O2  line  absorption  peaks. 


A  numerical  account  for  selected  frequencies  in  the  45-  to  75  GHz  range  is  given  in  Table  11.  Maximum 
shielding  (A2  >  150  dB)  of  satellite-based  systems  against  interference  from  ground  levels  can  be  acnieved 
between  58.9  and  61.2  GHz.  The  curve  of  associated  phase  dispersion  (Eq.  21)  is  not  shown.  Frequencies 
around  57  and  63  GHz  yield  max/min  values  of  .im  =  +  14  radians,  while  the  nondispersive  delay  of  the  air 
mass  is 


tin)  ’  7.64  +  0.21V  ns 

Both  Az  and  l4>z|  increase  with  decreasing  average  temperature  of  the  air  mass  (see  x  and  y,  Eq.  23) 
as  is.  the  case  for  different  climates. 

As  a  prospective  system  operates  at  increased  altitudes,  the  shielding  effect  of  the  unstructured  O2 
band  shape  breaks  down  and  about  30  communication  channels  open  up  between  the  lines  (see  Fig.  2d,  e). 

The  special  transfer  characteristics  of  these  channels  can  be  evaluated  by  applying  the  formulations  given 
in  Sections  2  to  4. 

Slanted  radio  paths  with  elevation  angles,  •p  >_  6°  (from  horizontal),  follow  a  secant  law  in  clear  air, 
i.e. , 


A^,  =  Az/sin'f  dB  (30) 

Refinements  at  low  elevation  angles  consider  refractive  and  dispersive  ray  bending  [H0PP0NEN,  1980], 


6.  PULSE  PROPAGATION 

Broadband  signal  transmission  at  carrier  frequencies,  fc  =  45  to  75  GHz,  through  the  clear  atmosphere  has 
to  reckon  with  the  radio  path  transfer  function  'll  (Eq.  22).  Digital  communication  systems  with  high  data 
rates  (0.5  to  3  Gbit/s)  are  possible,  employing  pulsed  signals  with  durations  xs  in  the  range  2  to  0.3  ns 
and  requiring  transmission  bandwidths,  EW,  between  2  and  >  10  GHz.  The  chances  for  pulse  shape  distortion 
increase  when  maximum  path  lengths  L  are  considered,  based  upon  path  attenuation  A.  For  examole,  a 
terrestrial  (h=0  km)  system  configuration  designed  to  tolerate  A  -  25  dB  could,  in  principle,  operate  at 
50  GHz  over  50  km,  which  reduces  to  1.6  km  at  60  GHz. 

Calculation  of  pulse  distortion  is  accomplished  as  a  function  of  pulse  shape,  pulse  duration,  and  path 
length  by  utilizing.  In  a  mini -computer,  the  method  of  Fast  Fourier  Transform  (FFT).  The  formulation 
starts  with  an  initial  signal  waveform 

eQ(t)  *  yr  /  EQ(f)exp(j2hft)df  (31) 

*  -OO 

where  e(t)  and  E0(f)  constitute  a  voltage  signal  Fourier  transform  pair.  The  signal  is  assumed  to  propagate 
as  a  plane  wave  and  Is  received  at  path  length  L  as 

e(t.L)  -  ^7  E0(f)TL(f)exp(j2sft)df  (32) 

The  nondispersive  propagation  time  (Eqs.  4  and  6) 

t°  *  3.336(106  +  NQ)l  ps  (33) 

is  factored  out  of  TL  and  numerical  evaluations  of  Eq.  (32)  are  performed  with  the  FFT  method  [ALLEN  and 
LIEBE,  1982],  Fields  of  individual  spectral  components  propagate  with  their  own  attenuation  and  velocity, 
and  their  superposit  on  at  the  receiver  may  lead  to  considerable  pulse  shape  distortions.  Higher  frequency 
components  of  a  pulse  signal  are  sped  up  with  respect  to  the  fundamental  for  negative  gradients  of  the 
phase  rate  dd/df  and  vice  versa.  The  received  pulse  might  experience  shortening  or  lengthening  of  its 
original  duration;  also,  it  can  break  up  into  several  parts  bearing  no  resemblance  to  the  original  signal, 
which  leads  to  intolerable  bit-error  rates. 

Propagation  properties  of  sea  level  air  are  enumerated  in  Table  III.  A  more  detailed  analysis,  covering 
the  altitude  range  h  *  0  to  30  km,  was  reported  by  LIEBE  and  HOPPONEN  (1977). 

Table  III.  Propagation  Properties  of  Moist  Sea  Level  Air 
(0  km,  15°C,  50*  RH) 


f 

GHz 

45.0 

50.0 

53.9 

55.8 

56.3  60.4 

63.0 

63.9 

65.8 

70.0 

75.0 

t 

dB/km 

0.18 

0.39 

2.11 

6.37 

7.85  15.05 

10.70 

7.42 

2.58 

0.59 

0.37 

3 

rad/ km 

0.22 

0.40 

0.81 

1.01 

0.99  -0.17 

-1.17 

-1 .25 

-1.03 

-0.58 

-0.40 

d./df 

dB/km  GHz 

0.019 

0.093 

1.31 

3.02 

3.10  0.00 

-3.90 

-3.61 

-1 .46 

-0.098 

-0.017 

♦MAX 

-MAX 

dd/df 

rad/km  GHz 

0.023 

0.060 

0.15 

0.00 

-0.03  -0.41 

-0.22 

0.00 

0.16 

0.057 

0.024 

-mAx 

tmSx 

Sea  level 

air  presents 

the  most 

favorable  conditions 

for  short-pulse  propagation.  Frequency  gradients  of 

individual  Oj  lines  at  higher  altitudes  can  increase  by  as  much  as  a  factor  of  100  on  a  MHz  frequency  scale. 
Examples  are  given  in  Fig.  7  for  three  types  of  pulse  signals  operated  over  a  link  with  the  design  parameters 
h  »  0  km,  clear  air,  L  *  35  km,  fc  *  *0  GHz,  r-  ■  1  ns.  The  signals  considered  are: 

(a)  Gaussian  shape  requiring  minimum  BW  («*  2  GHz). 
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(b)  Rectangular  shape  requiring  maximum  BW  («  +  10  GHz). 

(c)  Binary  phase  Shift  key  (BPSK)  modulation. 

Examples  of  more  serious  pulse  distortions  and  details  of  the  calculation  procedure  have  been  discussed 
by  ALLEN  and  LIE8E  P982].  It  is  convenient  to  include  in  Eq.  (32)  the  complex  frequency  response  of  a 
B'.J- 1  tmi  ting  device  (e.g.,  receiver  bandwidth,  smoothing  filter,  etc.),  or  a  UW-1  inearizinq  device  (e.g., 
adaptive  pnase  correction)  since  average  02  properties  [a(p,Tj,  .’.(p.T)]  for  a  given  application  might  De 
constant. 


7.  CONCLUSIONS 

Moist  air  has  been  characterized  as  an  atmospheric  propagation  medium.  The  basic  physics  of  molecular 
absorption  in  a  radio  path  have  been  cast  into  a  model  with  optimum  computer  run  time,  but  without  undue 
approximations.  The  paper  provides  graphical  and  numerical  examples  on  single  frequency  and  broadband 
behaviors  of  typical  radio  paths  operating  in  the  45  to  75  GHz  frequency  range.  The  routines  are  operated, 
even  by  nonexpert  users,  simply  by  calling  the  program  name,  choosing  from  a  battery  of  10  model  atmospheres, 
and  typing  in  the  interactively  requested  input  parameters  —  meteorological  data  are  converted  into 
propagation  data. 

The  program  is  available  in  FORTRAN  IV  upon.request,  together  with  print-outs  at  selected  freouencies 
[f  »  45,  46,  47,...  GHz;  f  »  Vg;  and  f  =  (v0  +  uq)/2]  over  48  height  levels  h  =  0  to  30  km  for  a  model 
atmosphere  (arctic-to-tropical  > .  Displayed  are  values  of  o,  da/df,  3,  d3/df,  A2,  :z,  V,  VJ,  and  t? 
for  relative  humidities,  RH  =  0  to  1003,  specified  over  h=0  to  8  km.  Cloud  and  statistical  rain  atten¬ 
uation  models  are  also  available.  Shortcomings  of  the  oresent  model  lie  in  the  omitted  spectroscom;  data 
base  for  trace  gases  and  in  the  empirical  nature  of  the  water  vaoor  contribution  r:  (161. 

Detrimental  effects  such  as  fluctuations  in  amplitude  [MEDEIROS-FILHO  et  al.,  1981],  phase,  and  group 
velocities,  and  the  direction  of  the  radio  wave,  may  be  predicted  from  spatial  and  temporal  variations  of 
p-T-e-w  distributions  along  a  radio  path.  A  priori  assessments  can  be  made  of  rare  events  for  reliable 
system  operations  and  of  possible  bandwidth  limitations. 
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Table  I.  Selected  Specific  Attenuation  ( i)  and  Phase  (s)  Rates  at  Altitudes 
h=0  to  30  km  for  U.S.  Std.  Atm.  [NOAA,  1976]. 


GHZ 

vo 

60.306 

5- 

60.400 

60.435 

7+ 

60.800 

61.151 

9+ 

61 .500 

61.800 

11  + 

62.100 

a 

15.1 

15.1 

15.1 

dB/km 

rad/km 

15.1 

14.9 

14.6 

14.1 

13.6 

8 

-  0.119 

-  0.145 

-  0.172 

-  0.318 

-  0.464 

-  0.599 

-  0.736 

-  0.364 

a 

13.7 

13.8 

13.8 

13.6 

13.3 

13.1 

12.7 

12.1 

8 

-  0.126 

-  0.152 

-  0.179 

-  0.321 

-  0.454 

-  0.574 

-  0.699 

-  0.821 

a 

11.3 

11.3 

11.3 

10.7 

10.3 

9.93 

9.72 

9.29 

6 

-  0.124 

-  0.157 

-  0.191 

-  0.329 

-  0.420 

-  0.497 

-  0.583 

-  0.633 

a 

9.51 

9.58 

9.48 

7.74 

7.60 

6.83 

7.05 

6.74 

6 

-  0.085 

-  0.145 

-  0.212 

-  0.339 

-  0.372 

-  0.409 

-  0.444 

-  0.499 

a 

6.56 

6.93 

6.70 

2.86 

4.63 

2.51 

4.14 

2.67 

8 

0.067 

-  0.067 

-  0.213 

-  0.222 

-  0.200 

-  0.213 

-  0.207 

-  0.203 

a 

4.52 

4.79 

4.89 

0.71 

3.71 

0.63 

3.23 

0.69 

8 

0.140 

-  0.011 

-  0.196 

-  0.112 

-  0.093 

-  0.100 

-  0.092 

-  0.030 

a 

3.20 

1.65 

3.78 

0.14 

3.39 

0.10 

2.94 

0.21 

8 

0.083 

0.006 

-  0.105 

-  0.043 

-  0.034 

-  0.036 

-  0.033 

-  0.039 

a 

2.91 

0.57 

3.38 

0.03 

3.28 

0.03 

2.86 

0.03 

6 

0.047 

0.004 

-  0.056 

-  0.022 

-  0.017 

-  0.019 

-  0.017 

-  0.014 

Attenuation  a  in  decibels  per  kilometer  and  refractive  dispersion  0  in  parts  per  million 
for  dry  (RH  *  0T)  to  saturated  (RH  *  100%)  air  at  sea  level  (h  -  0  km)  for  the  U.S. 
Standard  Atmosphere  [flOAA,  1976]  over  the  frequency  range  f  *  l  to  1000  GHz  [IIE8E,  1981] 


Table  II.  evaluation  of  Zenith  Attenuation  A  for  U.S.  Std.  Atm.  76  (Water  Vapor  V=15  mm) 
Between  45-  and  75  GHz,  Including  Mesospheric  Zeeman  effects 


(Figs.  2.3) 
h=0  (  30  km 


(Fig.  5) 
h=0  to  30  km 


A^ (0.3)  A 1 (0.6) 

(Fig.  4) 
h=30  to  100  km 


(Fig.  6) 
h=0  to  100  km 


Figure  3. 
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Attenuation  t  in  decibels  per  kilometer  and  refractive  dispersion  D  In  parts  per  million 
in  the  oxynen  microwave  band  at  altitudes  h  •  0-30  km  (U.S.  Std.  Atm.). 

a)  f  *  53.8  to  56.2  GHz  ,  b)  f  *  60.2  to  60.6  GHz  . 


Attenuation  a,  dB/km  Attenuation 


Fioure  4. 
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Zeeman  attenuation  patterns  of  the  oxygen  microwave  lines  K  -  19’  (see  Fig.  3a)  and  7 
(see  Fin.  3b)  for  altitudes  h  *  30-100  km  (U.S.  Std.  Atm.).  Each  frame  displays  -t(i,), 
o*(n,),  and  n*(.j,)  patterns  for  the  mannetic  field  strenaths  H  =  0.3  (left  hand)  and 
H  *  0.6  G  (right  hano).  The  patterns  are  symmetric  with  respect  to  the  center  axis  (inter- 
chanoe  and  c~).  The  frequency  deviation  A  *  (f  -  v  )  fs  between  t40  MHz  for  h  =  30  and 

t2  MHz  for  h  1  60-100  km.  The  maximum  attenuation  rate  n  is  that  of  the  isolated  unsplit 
line  (H  -  0  and  f  «  v  ). 


100 


Figure  5. 


Figure  6 


Zenith  attenuation  A  in  decibels  of  the  U.S.  Std.  Atm.  for  a  cl oud-free, dry  (V  5  1.5  mm) 
moist  (15.3  mm) ,  and  saturated  {30.5  mm,  PH  -  100"  for  h  *  0  to  8  km)  air  mass  over  the 
'requency  ranoe  f  *  10  to  350  GHz, 


retails  of  zenith  attenuation  A  ( V  =  15.  m,  see  Fig.  5)  over  the  frequency  range  f  *  53 
67  GHz  exhibiting  60  GHz  oxygen  absorption  with  band  and  complete  line  structure.  A 
m-type  Zeeman  pattern  (t,)  (see  Fig.  4)  and  H  ■  0.3  Gauss  were  assumed  for  the  height 
range  h  *  30  to  100  Vm.  1 
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Figure  7.  Relative  pulse  shape  distortion  over  a  terrestrial  path  of  35  km  lei 
frequency  of  50  GHz  using  Gaussian,  rectangular,  and  UPSK  1-ns-wide 
(dashed  line). 
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